Small-Caliber  Projectile  Target  Impact 
Angle  Determined  from  Close 
Proximity  Radiographs 


by  Tyler  E.  Ehlers,  Bernard  J.  Guidos,  David  W.  Webb 


ARL-TR-3943 


October  2006 


Approved  for  public  release;  distribution  is  unlimited. 


NOTICES 

Disclaimers 

The  findings  in  this  report  are  not  to  be  construed  as  an  official  Department  of  the  Army 
position,  unless  so  designated  by  other  authorized  documents. 


Citation  of  manufacturers’  or  trade  names  does  not  constitute  an  official  endorsement  or 
approval  of  the  use  thereof. 

DESTRUCTION  NOTICE — Eor  classified  documents,  follow  the  procedures  in  DoD  5220.22- 
M,  National  Industrial  Security  Program  Operating  Manual,  Chapter  5,  Section  7,  or  DoD 
5200. 1-R,  Information  Security  Program  Regulation,  C6.7.  Eor  unclassified,  limited  documents, 
destroy  by  any  method  that  will  prevent  disclosure  of  contents  or  reconstruction  of  the  document. 


Army  Research  Laboratory 

Aberdeen  Proving  Ground,  MD  21005-5069 


Approved  for  public  release;  distribution  unlimited. 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
0MB  No.  0704-0188 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing 
data  sources,  gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  information.  Send  comments  regarding  this  burden  estimate 
or  any  other  aspect  of  this  collection  of  information,  including  suggestions  for  reducing  the  burden,  to  Department  of  Defense,  Washington  Headquarters  Services, 
Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302.  Respondents  should  be  aware 
that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a  collection  of  information  if  it  does  not  display  a 
currently  valid  0MB  control  number. 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 


1.  REPORT  DATE  (DD-MM-YYYY) 

October  2006 


2.  REPORT  TYPE 


Final 


3.  DATES  COVERED  (From  -  To) 

Jan  2005  -  Jun  2006 


4.  TITLE  AND  SUBTITLE 


Small-Caliber  Projectile  Target  Impact  Angle  Determined  from  Close  Proximity 
Radiographs 


5a.  CONTRACT  NUMBER 

NA 


5b.  GRANT  NUMBER 


5c.  PROGRAM  ELEMENT  NUMBER 


6.  AUTHOR(S) 

Tyler  E.  Ehlers,  Bernard  J.  Guidos,  David  W.  Webb 


5d.  PROJECT  NUMBER 

62.AH80/AH84 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

U.S.  Army  Research  Laboratory 
ATTN:  AMSRD-ARL-WM-TC 
Aberdeen  Proving  Ground,  MD  21005-5069 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

ARL-TR-3943 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

U.S.  Army  Research  Laboratory 
Aberdeen  Proving  Ground,  MD  21005-5069 


10.  SPONSOR/MONITOR'S  ACRONYM(S) 


11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  is  unlimited. 


13.  SUPPLEMENTARY  NOTES 


14.  ABSTRACT 

Small-caliber  rounds  have  a  yaw  cycle  of  approximately  three  meters,  but  the  magnitude  and  initial  yaw  varies  from  shot  to  shot,  even  when 
fired  from  the  same  gun.  This  yaw  cycle  becomes  a  variable  that  researchers  must  consider  when  looking  at  the  lethality  of  an  individual 
shot  and  round  type,  especially  when  testing  at  close  ranges  where  the  yaw  cycle  is  high.  Four  orthogonal  radiograph  images  provided 
measured  pairs  of  pitch,  yaw  and  their  Z-locations,  a  curve  is  fitted  through  all  four  pair  of  points  using  a  fundamental  Aerodynamic 
equation.  The  fitted  curve  is  then  carried  out  a  known  distance  to  the  target  where  the  curve  provides  an  expected  pitch,  yaw,  and  total  angle 
of  attack  at  impact.  Upon  completion  of  a  significant  number  of  shots,  this  method  can  provide  an  accurate  estimate  of  measurement  error. 


15.  SUBJECT  TERMS 

small  caliber,  yaw  cycle,  fundamental  aerodynamic  equation,  error  analysis 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION 

OF 

ABSTRACT 

18.  NUMBER 

OF 

PAGES 

19a.  NAME  OF  RESPONSIBLE  PERSON 

Tyler  E.  Ehlers 

a.  REPORT 

b.  ABSTRACT 

c.  THIS  PAGE 

SAR 

51 

19b.  TELEPHONE  NUMBER  (Include  area  code) 

u 

u 

u 

410-278-1014 

Standard  Form  298  (Rev.  8/98) 

Prescribed  by  ANSI  Std.  Z39.18 


11 


Contents 


List  of  Figures  iv 

Acknowledgments  v 

1,  Introduction  1 

2,  Experiment  and  Analysis  1 

2.1  Range  Setup . 1 

2.2  Radiograph  Data  Reduction . 3 

2.3  Projectile  Motion . 7 

3,  Results  12 

4,  Error  Estimation  16 

4.1  Trajectory  Angle  Error . 16 

4.2  Target  Impact  Angle  Error . 22 

5,  Conclusions  28 

References  29 

Appendix  A:  Maximum  Likelihood  Estimation  of  the  Displacement  Standard  Deviation  31 

Appendix  B:  Expected  Value  of  the  Square  Root  of  an  Inverted  Gamma  Random  Variable33 

Appendix  C:  Maximum  Likelihood  Estimation  of  the  Three  Dimensional  Displacement 
Standard  Deviation  35 

Abbreviations  and  Symbols  37 


Distribution  List 


39 


List  of  Figures 


Figure  1.  Illustration  of  experimental  setup . 2 

Figure  2.  Example  of  radiograph  of  M855  projectile . 3 

Figure  3.  Projectile  geometry  template . 4 

Figure  4.  Magnification  factor  stand . 5 

Figure  5.  Digitizing  program  digitized  points  and  intermediate  calculated  values . 6 

Figure  6.  Digitizing  program  output . 7 

Figure  7.  Fast  and  slow  epicyclic  modal  arms . 8 

Figure  8.  6-DOF  simulation  of  M855  pitch  angle  vs.  yaw  angle . 9 

Figure  9:  Gamma  at  Impact  Prediction  Excel®  Spreadsheet . 11 

Eigure  10.  Measured  and  fitted  pitch  angle  vs.  yaw  angle . 12 

Eigure  1 1 .  Measured  and  fitted  pitch  and  yaw  angle  for  shot  21 . 13 

Eigure  12.  Measured  and  calculated  total  impact  angle  for  shot  21 . 14 

Eigure  13.  Distribution  of  measured  impact  yaw  occurrences,  M855,  101  shots . 15 

Eigure  14.  Comparison  of  relative  frequency  histograms  of  experimental  RMS  and  simulated 
RMS  under  assumption  of  constant  errors  in  the  horizontal  direction . 17 

Eigure  15.  Comparison  of  relative  frequency  histograms  of  experimental  RMS  and  simulated 
RMS  under  assumption  of  constant  errors  in  the  vertical  directions . 18 

Eigure  16.  Comparison  of  relative  frequency  histograms  of  experimental  RMS  and  simulated 
RMS  under  assumption  of  gamma-distributed  (a=.4027,  b=3.25)  errors  in  the  horizontal 
direction . 20 

Eigure  17.  Comparison  of  relative  frequency  histograms  of  experimental  RMS  and  simulated 
RMS  under  assumption  of  gamma-distributed  (a=.1684,  b=7.2)  errors  in  the  vertical 
direction . 20 

Eigure  18.  Estimated  standard  deviation  of  trajectory  as  a  function  of  reading  error . 21 

Eigure  19.  Least-squares  fit  RMS  distribution  from  Monte  Carlo  simulations . 23 

Eigure  20.  Comparison  of  relative  frequency  histograms  of  experimental  RMS  and  simulated 
RMS  under  assumption  of  constant  errors . 25 

Eigure  21 .  Comparison  of  relative  frequency  histograms  of  experimental  RMS  and  simulated 
RMS  under  assumption  of  gamma-distributed  (a=.0164,  b=12.25)  errors . 25 

Eigure  22.  Monte  Carlo  simulation  of  impact  yaw  error  distribution . 26 

Eigure  23.  Estimated  standard  deviation  of  total  impact  yaw  as  a  function  of  reading  error . 27 


IV 


Acknowledgments 


The  authors  gratefully  acknowledge  the  following  organizations  and  individuals. 

The  guidance  of  Lee  Magness  is  greatly  appreciated,  without  his  dedication  this  program  would 
not  have  been  possible.  The  authors  would  like  to  thank  the  Jim  Newill  for  funding  the  program. 
We  would  also  like  to  thank  Jason  Angel  for  initiating  the  Excel  Digitizing  Program,  Bryan 
Peterson  for  developing  the  templates  and  Rick  Summers  for  continuous  forethought  and 
direction  in  optimizing  the  Monte  Carlo  Visual  Basic®  program. 

Additionally,  the  authors  would  like  to  express  gratitude  to  the  government  technicians  Rick 
Kane,  Mark  Ward,  and  the  Dynamic  Science  Inc.  technicians  Maurice  Clark,  Mark  Hale,  and 
Ron  Berry  for  orchestrating  the  experiments.  It  is  with  great  appreciation,  that  we  thank  Greg 
Watt  for  digitizing  every  shot.  Without  his  dedication,  this  project  would  not  have  been 
successful. 


V 


Intentionally  Left  Blank. 


VI 


1.  Introduction 


For  small-caliber  projectiles  sueh  as  the  5.56-millimeter  (mm)  M855,  assessments  of 
terminal  ballistic  performance  versus  soft  targets  are  condueted  by  firing  a  single  round 
into  a  ballistic  gelatin  bloek.  The  gelatin  bloek  is  then  dissected  and  the  terminal  ballistic 
quantities  of  interest  are  obtained  by  careful  measurements  of  the  projectile’s  path  in  the 
target.  The  proeess  is  repeated  for  a  suflioient  number  of  shots,  and  the  results  are 
eompiled  to  form  a  terminal  ballistic  model. 

The  terminal  ballistic  quantities  of  interest  are  commonly  represented  by  mathematical 
functions  that  use  target  impact  angle  as  an  independent  variable.  The  target  impact 
angle  is  a  sealar  quantity  defined  here  as  the  angle  between  the  projectile  axis  and  a  line 
normal  to  the  target  impact  surface.  The  quality  of  the  terminal  ballistie  model  is 
influeneed  by  the  accuracy  of  the  measured  target  impaet  angle  from  eaeh  experiment. 
For  small-ealiber  projectiles,  accurate  determination  of  the  target  impact  angle  requires  a 
consideration  of  aero-ballistic  free-fiight  motion  model  parameters  and  awareness  of 
photographic  measurement  limitations. 

This  report  presents  the  methodology  reeently  developed  at  the  U.S.  Army  Research 
Laboratory  (ARL)  to  calculate  the  target  impact  yaw  from  data  measured  from  close 
proximity  radiographs  (i.e.,  “x-ray  images”).  The  experimental  setup  and  proeedures 
used  to  capture  the  radiograph  images  are  detailed  in  section  2.1.  Section  2.2  covers  the 
step-by-step  proeedures  required  to  reduce  the  radiograph  images  into  data  points.  Going 
further,  these  data  points  are  used  in  section  2.3  where  they  are  fit  to  the  analytical  aero- 
ballistic  free-flight  model.  The  aerodynamic  fit  is  then  extrapolated  from  the  last 
radiograph  to  the  target  to  aceurately  detail  the  impact  angle  of  the  projectile.  Seetion  3 
covers  the  results  of  testing  using  this  method,  and  how  the  data  sets  can  be  manipulated. 
Section  4  covers  a  procedure  for  caleulating  the  error  of  the  determined  target  impact 
yaw;  additionally  the  results  for  recent  small-caliber  firings  are  shown  and  the 
measurement  error  is  estimated. 


2.  Experiment  and  Analysis 


2,1  Range  Setup 

The  experimental  setup  followed  the  basic  principles  presented  in  Zook,  et  al.  (1992)  and 
is  illustrated  in  figure  1 .  The  gun  muzzle  was  positioned  between  5  and  7  meters  (m) 
from  the  gelatin  block  striking  surface  at  0°  quadrant  elevation.  Four  orthogonal  pairs  of 
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radiograph  stations  were  situated  to  capture  images  at  0.23-m  [9-inch  (in.)]  intervals  up- 
range  from  the  gel  block.  The  fourth,  most  down-range  radiograph  station  was  situated 
0.23  m  from  the  gel  block.  The  150  keV  x-ray  film  was  positioned  0.18  m  to  the  left  of 
the  line  of  fire.  The  150  keV  x-ray  heads  were  .50  m  to  the  right  of  the  line  of  fire. 
Therefore  the  distance  between  the  x-ray  generating  tubes  and  the  film  was  0.68  m. 


Vertical  Fiducial  Wire 
•4 - Gel  Block  fTarget) 


Figure  1.  Illustration  of  experimental  setup. 

The  side  tube  heads  were  manually  aligned  with  a  level  to  ensure  they  were  placed  on  a 
single  plane.  The  top  tube  heads  were  manually  aligned  with  a  plumb  bob  to  ensure  they 
were  placed  on  a  single  plane.  If  the  heads  are  aligned  carefully  during  this  process,  the 
tube  heads  should  be  within  2  mm  of  the  intended  plane. 

Prior  to  shooting,  two  24-gauge  wound-steel  fiducial  wires  were  suspended.  The 
radiograph  films  were  pressed  against  the  fiducial  wires,  which  serve  as  reference  points. 
These  reference  points  were  later  used  to  calculate  projectile  angular  orientation  and 
lateral  position  in  their  respective  orthogonal  planes.  The  side-film  fiducial  wire  was 
positioned  horizontally  parallel  to  the  boresighted  line  of  fire.  These  axially  aligned 
wires  were  manually  pulled  taut  (to  an  estimated  40-lb  tension)  and  clamped.  Because  the 
wire  length  was  only  about  one  meter,  it  is  assumed  that  no  significant  catenary  droop 
exists,  and  no  such  correction  was  applied  during  the  subsequent  analysis. 

For  each  radiograph  film,  a  wire  was  suspended  perpendicular  to  the  shot  line  fiducial 
wires  at  the  same  exact  downrange  location  as  one  of  the  radiograph  heads.  These 
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perpendicular  wires  detail  the  location  of  the  head  and  were  subsequently  used  as 
reference  points  when  calculating  the  projectiles  downrange  location.  The  system  of 
wires  was  kept  in  place  for  multiple  firings.  The  estimated  time  delay  for  each  station 
was  programmed  into  a  data  acquisition  system,  the  projectile  was  launched,  break 
screens  initiated  the  time  delay  generators,  and  radiographs  were  obtained  of  the  in-flight 
projectile  and  the  fiducial  wires. 

2.2  Radiograph  Data  Reduction 

A  total  of  four  1 1  in.  x  17  in.  radiograph  films  are  obtained,  each  having  two  images  of 
the  projectile,  as  shown  in  figure  2.  The  eight  projectile  images  are  measured  using  a 
light  table  and  computer  digitization  software.  These  measurements  are  then  processed 
to  yield  projectile  location  and  angular  orientation  relative  to  the  line  of  fire  at  the  four 
orthogonal  stations. 

Mi*m. 


Ssot 


Figure  2.  Example  of  radiograph  of  M855  projeetile. 

For  small-caliber  projectiles,  a  geometry  measurement  template  can  improve  the  speed 
and  accuracy  of  the  manual  digitization  process.  The  template  consists  of  a  transparency 
having  concentric  images  and  reference  points  of  the  projectile  at  various  incidence  in  the 
operator’s  field  of  view.  These  images  then  have  a  template  laid  over  them  to  accurately 
align  the  orientation  of  the  bullet  and  place  reference  points  on  the  film.  These  reference 
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points  are  located  on  the  corners  of  the  box  in  figure  3,  and  are  centered  over  the  center 
of  gravity  (c.g.)  of  the  bullet. 


— 

Figure  3.  Projectile  geometry  template. 


The  specific  digitizing  process  is  as  follows.  A  single  piece  of  film  is  placed  on  a  light 
table  with  the  bullet  flying  right  to  left.  Using  a  light  table  and  Visual  Basic®  digitizing 
program,  the  user  takes  measurements  of  the  locations  of  the  fiducial  wires,  a  fixed  point 
on  the  bullet,  and  the  four  reference  points,  which  the  user  must  digitize  in  the  following 
order.  The  first  mark  is  the  crossing  of  the  fiducial  wires  on  the  right-hand  side  of  the 
film.  The  second  mark  is  located  where  the  user  can  accurately  place  the  cursor  on  the 
extreme  left  side  of  the  fiducial  wire.  The  third  mark  is  the  tip  of  the  bullet  of  the  first 
flash  which  is  the  right  image.  The  fourth  through  seventh  marks,  marks  made  using  the 
template,  are  the  top  left  reference  mark,  followed  by  the  bottom  left  reference  mark, 
followed  by  the  top  right  mark,  concluding  with  the  bottom  right.  The  second  image  is 
digitized  in  the  same  manner  as  the  first  bullet  image.  The  entire  process  is  repeated  for 
each  projectile  image.  The  Visual  Basic®  program  outputs  the  measurements  to  an 
Excel®  spreadsheet. 

These  raw  measurements  are  input  into  the  digitizing  Excel®  program,  shown  in  figures 
5  and  6,  which  calculates  the  projectile  angular  orientation.  The  angle  in  the  vertical 
plane  (i.e.,  nose  pointed  up  or  down)  is  the  pitch  angle  a.  The  angle  in  the  horizontal 
plane  (i.e.,  nose  pointed  right  or  left)  is  the  yaw  angle  ji.  The  total  angle  of  attack  is  y  • 

The  program  calculates  the  magnification  factor,  also  referred  to  as  K-factor,  so  the 
actual  location  of  the  bullet  when  the  image  was  captured  can  be  found.  This  program 
also  calculates  the  trajectory  angle,  i.e.,  the  angle  that  the  projectile  is  traveling  in  relation 
to  the  axial  fiducial  wires.  The  trajectory  angle  component  (i.e.,  the  angle  between  the 
projectile  velocity  vector  and  the  fiducial  wire)  in  the  vertical  plane  is  denoted  rj^  and  the 
trajectory  angle  component  in  the  horizontal  plane  is  denoted  rj p  . 

The  fiducial  wires  are  digitized  and  used  as  the  reference  plane.  Because  the  fiducial 
wires  set  the  reference  plane,  the  radiograph  can  be  positioned  on  the  light  table  at  any 
angle.  The  projectile’s  reference  points  are  digitized  to  locate  the  centerline  of  the  bullet 
and  coordinates  of  the  image.  Eor  both  orthogonal  views,  the  angle  of  the  bullet  relative 
to  the  fiducial  wire  is  found  for  each  image  using  the  bullet  centerline  and  the  fiducial 
wire.  The  angle  of  the  bullet  relative  to  the  fiducial  wire  is  referred  to  here  as  the 
apparent  angle.  The  apparent  projectile  pitch  angle,  aA,  is  determined  from  the  side  fdm; 
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the  apparent  yaw  angle,  Pa,  is  determined  from  the  bottom  film.  The  sign  conventions 
used  for  pitch  and  yaw  angles  are  positive  nose  up  and  to  the  right. 

To  calculate  the  velocity  and  the  actual  location  of  the  projectile,  the  magnification  factor 
must  be  determined.  The  derivation  of  the  magnification  factor  is  given  in  Zook,  et  al. 
(1992,  section  5.1);  the  pertinent  equation  is  restated  here: 

MagnificationFactor  =  ^ HFoid hf  -  d -  {d -  d ao  -  d 

{.d  HFod  HF  )  “  i.^AO  ~  dffO  )(.^A  ~  d  ff  ) 


Where: 


duF  is  the  distance  from  the  head  to  the  film, 

dnFO  is  the  distance  from  the  head  to  the  film  in  the  orthogonal  view, 

dFF  is  the  distance  from  the  fiducial  wire  to  the  film, 

dFFO  is  the  distance  from  the  fiducial  wire  to  the  film  in  the  orthogonal  view, 

Ia  is  the  apparent  location, 

Iao  is  the  apparent  location  in  the  orthogonal  view. 


In  order  to  model  the  projectile  angular  motion,  the  pitch  and  yaw  angles  relative  to  the 
velocity  vector,  denoted  a  and  P,  respectively,  are  needed.  They  are  obtained  by 
correcting  ou  and  Pa  by  the  trajectory  angles  in  both  orthogonal  planes,  rj^  and//^  , 

respectively.  The  values  of  and//^  are  calculated  by  taking  a  linear  least  squares  fit  of 
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the  projectile  c.g.  lateral  location  at  each  of  the  four  radiograph  stations  as  a  function  of 
range. 

The  values  of  and  7^  are  assumed  constant  over  the  flight  interval  of  interest.  Linear 

aerodynamic  theory  demonstrates  that  the  projectile  trajectory  experiences  an  epicyclical 
(or  corkscrew-like)  motion,  referred  to  as  the  swerving  motion  (e.g.,  McCoy,  1998).  The 
trajectory  angle  components,  7^  and 7^  ,  vary  in  a  fluctuating  manner  with  respect  to 

range.  The  magnitude  of  the  fluctuation  in  each  plane  is  easily  calculated  using  linear 
aerodynamic  theory,  and,  for  the  parameters  and  conditions  of  this  study,  is  insignificant. 

For  small  angles  of  attack,  the  total  angle  of  attack  can  be  approximated  as; 

r  =  ■  (1) 


When  testing  small-caliber  ammunition,  as  a  and  fi  are  often  greater  than  2  degrees  (deg), 
the  above  approximation  loses  accuracy.  In  the  present  study,  the  following  exact 
relationship  is  used; 
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tan 


tan^  a  +  tan^  /? 


(2) 
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Figure  5.  Digitizing  program  digitized  points  and  intennediate  caleulated  values. 
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Figure  6.  Digitizing  program  output. 

At  this  point  the  data  has  been  redueed  to  a,  ji,  and/  in  four  different  locations,  those 
four  locations  in  three  dimensional  space,  and  the  velocity  of  the  projectile. 


2,3  Projectile  Motion 

The  target  impact  angle  is  determined  from  a  mathematical  fit  of  the  discrete  data 
measured  in  the  radiographs  to  the  linear  aerodynamic  free-fiight  equations  of  motion  as 
outlined  by  McCoy  (1998).  The  complex  angle  of  attack  relative  to  the  flight  path  is 
defined  as: 


=  a  +  ip ,  (3) 

with  the  pitch  angle  a  being  positive  for  nose  up  and  the  yaw  angle  being  positive  for 
nose  right  from  the  gunner’s  perspective. 

The  projectile  is  assumed  to  be  symmetric,  spinning,  and  in  flat  fire.  The  modeled 
trajectory  length  is  assumed  to  be  of  low  angle  of  attack  and  short  enough  to  neglect 
damping.  Under  these  assumptions,  the  behavior  of  the  complex  angle  of  attack  is  given 
in  “modal  form”  by: 


i <Pfo  J ‘l>so  p  ‘t>'s  (^-^0 ) 


(4) 


in  which  x  is  the  distance  along  the  trajectory,  is  the  reference  location,  and  the  modal 
parameters  are  defined  as  follows: 
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Kp  =  magnitude  of  the  fast  epieyelie  modal  arm,  real  and  positive, 

=  magnitude  of  the  slow  epieyolie  modal  arm,  real  and  positive, 

(/)p^  =  referenee  phase  angle  of  the  fast  epieyelie  modal  arm,  evaluated  at  x=  xo  and  real, 
=  referenee  phase  angle  of  the  slow  epieyelie  modal  arm,  evaluated  at  x=  xo  and  real, 
(p'p  =  turning  rate  of  the  fast  epieyelie  modal  arm,  real, 

(pg  =  turning  rate  of  the  slow  epieyelie  modal  arm,  real. 

As  viewed  from  an  onlooker  standing  at  the  muzzle,  the  angular  motion  looks  like  the 
nose  of  the  projeetile  was  etehing  a  eircular  motion;  the  slow  arm,  while  simultaneously 
etehing  another  eireular  motion  about  that  eireular  motion;  the  fast  arm.  Figure  7  depiets 
this  motion.  As  the  bullet  travels  away,  both  arms  turn  in  a  eloekwise  motion  as  the 
projeetile  spins  in  a  eloekwise  direction. 


Figure  8  depicts  the  motion  of  the  tip  throughout  a  half  of  the  slow  arm  cycle. 
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Alpha  vs  Beta 


Figure  8.  6-DOF  simulation  of  M855  pitch  angle  vs.  yaw  angle. 

The  pitch  and  yaw  angles  can  be  expressed  separately  as: 

a  =  K,  cos[^/i^„  +  (x  -  x^)]  +  cos[^5„  +  ^5  “  ^0)]  (5) 

15  =  Kp  sin[^^„  +  ^;  (x  -  Xg)]  +  sin[4„  +  (x  -  x^)] .  (6) 

For  a  single  shot,  digitization  and  reduction  of  the  radiographs  produces  the  following 
known  quantities:  x, ,  a,. ,  p. ,  for  i  =  in  which  the  subscript  i  is  the  radiograph 

station  out  of  n  total  stations.  The  reference  location  Xg  is  taken  as  the  midpoint  of  the 

trajectory  encompassed  by  the  radiographs.  The  six  unknown  variables  are  the  modal 
parameters:  Kp  ,Kg ,  (p'p  ,  and  (f>'^  . 

The  turning  rates,  </)p  and  (p'^  ,  are  a  function  of  the  projectile’s  physical  properties,  spin 

rate,  velocity,  atmospheric  conditions,  and  the  aerodynamic  pitching  moment  coefficient 
( ^Ma  )•  The  pitching  moment  coefficient  is  a  measure  of  the  aerodynamic  moment  in  the 
plane  of  the  instantaneous  angle  of  attack.  For  the  M855  class  of  ammunition,  is 
positive,  indicating  that  spin  is  required  to  stabilize  the  bullet  in  flight.  The  values  of  p'p 
and  (p'g  (and  thus  )  have  traditionally  been  determined  from  spark  range 
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experimentation,  whieh  uses  spark  shadowgraph  imagery  to  provide  diserete  motion  data 
that  can  be  numerically  modeled  using  linear  aerodynamic  theory  or  6-degrees-of- 
freedom  equations  of  motion.  The  values  of  (l)'p  and  cannot  be  accurately  obtained 

from  the  four  radiograph  stations  used  in  the  present  study  because  the  trajectory  length  is 
too  short  and  the  number  of  stations  is  too  small.  Instead,  it  is  assumed  that  the  turning 
rates  have  already  been  accurately  determined.  By  prescribing  the  turning  rates,  the 
number  of  unknown  variables  for  each  shot  is  reduced  to  four:  Kp  ,  ,  (/)p^  ,  and  (j)^^ . 

When  fitting  equations  5  and  6  to  the  experimental  data,  it  is  necessary  to  specify  realistic 
initial  values  for  the  four  parameters  to  be  fitted.  In  this  study,  the  fitting  procedure  was 
incorporated  into  a  spreadsheet  that  included  a  graph  of  a  versus  (3,  to  facilitate  the 
specification  of  initial  values  and  verify  the  integrity  of  the  fitted  parameters.  The  four 
radiograph  stations  encompassed  approximately  14  of  a  fast  mode  cycle  length,  so  the 
slow  mode  magnitude  and  orientation  were  easy  to  visualize  as  the  center  of  a  nearly 
circular  pattern  comprised  of  the  four  data  points.  As  such,  it  was  convenient  to  specify 
the  center  of  this  circular  pattern  using  reference  pitch  and  yaw  angles,  denoted  and 

,  rather  than  specifying  and  (j)^^ .  The  trigonometric  relationships  are: 

Ks  =  (7) 

A.  =  (8) 

An  iterative  technique  is  used  to  fit  equations  5  and  6  to  the  experimental  data.  The 
technique  determines  the  values  of  the  four  unknown  variables  by  minimizing  the  root- 
mean-square  (RMS)  error  of  the  fit.  The  RMS  error  of  the  fit  is  defined  as: 


“^■5= 

V  i=\  i=\ 

in  which 

=ai-a(x,) 

=  Pi- 


(9) 


(10) 

(11) 


In  these  above  expressions,  «,  and  Pi  are  the  values  of  pitch  angle  and  yaw  angle 
determined  from  the  radiograph,  and  a(xi )  and  P(x , )  are  the  values  of  pitch  angle  and 
yaw  angle  determined  from  the  fitted  parameters  and  evaluated  at  x  / ,  the  location  of  the 
radiograph. 
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After  the  four  fitted  parameters  are  determined,  the  modeled  pitch  and  yaw  angles  at  the 
target,  ar  and  are  calculated  by  substituting  the  values  back  into  equations  5  and  6 
and  solving  at  xt.  The  impact  angle  of  the  projectile  at  the  target,  /t,  is  found  from  the 
equation: 


/j,  =  tan  '  -yjtan^  +  tan^  Pj.  .  (12) 

In  practice,  the  user  imports  the  calculated  in-flight  values  from  the  digitizing  program 
into  the  Gamma  at  Impact  Prediction  Excel®  Spreadsheet  shown  in  figure  9.  The 
Gamma  at  Impact  Spreadsheet  uses  the  linear  aerodynamics  model,  equations  5  and  6,  to 
calculate  ar  and  Pt  which  are  used  to  calculate  the  yr.  The  six  modal  parameters  are 
compiled  in  the  right-most  column. 
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Figure  9:  Gamma  at  Impact  Prediction  Excel®  Spreadsheet. 

To  run  the  iterative  error  reducing  program: 

1)  Click  Tools  >Add-Ins  >Solver, 


Calculated 
RMS  Error 


The  Solver  will  then  be  available  to  open. 
2)  Click  Tools  >Solver, 


The  Solver  will  issue  a  menu  box. 


3)  Select  the  Target  Cell  equal  to  the  calculated  RMS  Error  (pointed  out  in 
figure  9). 

4)  Set  the  Equal  To  box  to  Minimize, 

Set  the  By  Changing  Cells  to  contain  the  cells  containing  the  values  of  Kfast,  phif_0, 
Alpha_center,  and  Beta  center. 
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3.  Results 


The  fitted  results  shown  in  the  spreadsheet  of  figure  9  were  obtained  by  preseribing  the 
values  of  the  fast  arm  turning  rate,  (/)'p  ,  to  be  21 1  degrees  per  meter  (deg/m),  and  the 

slow  arm  turning  rate,  (j)'^ ,  to  be  30  deg/m.  These  values  were  obtained  from  the  spark 

range  results  reported  by  McCoy  (1985)  and  verified  by  more  recent  unpublished  in- 
house  spark  range  firings. 

Figure  10  shows  a  typical  result  of  the  measured  pitch  angle  versus  yaw  angle,  comparing 
the  actual  data  to  the  modeled  motion  data.  The  four  radiographs  encompass 
approximately  14  of  a  fast  mode  cycle  length.  The  extrapolated  target  impact  values  are 
also  plotted,  and  located  approximately  14  of  a  fast  mode  cycle  downrange  from  the  first 
radiograph  station.  Figure  1 1  shows  the  same  data  plotted  versus  range,  illustrating  the 
sinusoidal-like  motion  of  the  projectile  in  the  pitch  and  yaw  planes. 


M855  Shot#  21 


P  (deg) 


Figure  10.  Measured  and  fitted  pitch  angle  vs.  yaw  angle. 
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M855  Shot#  21 


Range  (m) 


Figure  1 1 .  Measured  and  fitted  piteh  and  yaw  angle  for  shot  21 . 

Figure  12  shows  the  measured  total  angle,  eomparing  the  aetual  data  to  the  modeled 
motion  data.  The  extrapolated  value  at  the  target,  yj ,  is  apparent  in  this  view.  This 
figure  illustrates  that  using  a  linear  fit  of  total  angle  obtained  from  the  radiograph  data 
has  the  potential  to  be  highly  inaeeurate.  Aeeurate  determination  of  target  impaet  angle 
using  a  linear  fit  of  total  angle  would  require  the  stations  to  be  elustered  extremely  elose 
to  the  target.  The  approach  may  be  unfeasible  considering  that  the  current  positioning  of 
radiograph  equipment  already  uses  the  available  space. 
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M855  Shot#  21 


Range  (m) 


Figure  12.  Measured  and  caleulated  total  impact  angle  for  shot  21. 

Figure  13  shows  the  target  impact  yaw  for  101  shots  of  the  M855  projectile  from  the 
same  barrel.  The  result  is  presented  as  a  distribution  from  0  to  5  deg  in  0.5  deg 
increments.  For  the  relative  muzzle  and  target  positions  used  for  these  firings,  the  most 
common  range  of  total  impact  angle  values  was  1.0  to  1.5  deg  (25  occurrences).  The 
second  most  common  range  was  0.5  to  1.0  deg  (20  occurrences).  Considering  that  the 
range  0  to  0.5  deg  had  1 1  occurrences,  56  of  101  shots  had  impact  angles  less  than  2  deg. 
The  largest  total  impact  angle  observed  during  the  firings  was  between  5.0  and  5.5  deg, 
with  one  occurrence.  This  series  of  data  shows  how  a  data  set  can  be  manipulated  by 
repositioning  the  gun  barrel  to  capture  the  range  of  yaws  that  are  commonly  found  in 
CQB  environments. 
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M855  Yaw  at  Impact  vs.  Frequency  Observed 
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Figure  13.  Distribution  of  measured  impaet  yaw  oecurrenees,  M855,  101  shots. 

By  monitoring  the  development  of  the  impact  angle  distribution  during  a  firing  program, 
the  distribution  of  target  impact  angle  occurrences  could  be  tailored  for  specific  purposes. 
Such  tailoring  might  be  desirable;  for  example,  to  obtain  more  occurrences  of  yaw  levels 
greater  than  2  deg  in  order  to  more  accurately  characterize  the  terminal  ballistics  in  that 
regime.  The  tailoring  could  be  achieved  by  adjusting  the  distance  between  the  muzzle 
and  the  target,  using  the  total  yaw  behavior  illustrated  in  figure  12  as  observed  for  a  large 
sample  of  shots.  Concurrently,  usable  estimates  of  desired  muzzle-to-target  separation 
distance  that  minimizes  or  maximizes  the  target  impact  angle  could  be  achieved  using  the 
turning  rates  of  the  projectiles. 

As  an  illustrative  example,  suppose  the  experimentalist  completed  the  fit  from  shot  21, 
figures  10,  11,  and  12,  and  wanted  to  capture  the  minimum  yaw  on  shot  22.  Figures  1 1 
or  12  indicate  that  by  repositioning  the  muzzle  0.68  m  closer  to  the  target  (shifting  the  fit 
to  the  right  on  the  graph)  for  shot  22,  the  projectile  will  be  more  likely  to  impact  the 
target  during  a  portion  of  its  yaw  cycle  where  minimum  values  occur.  On  the  other  hand, 
suppose  the  experimentalist  wanted  to  capture  the  maximum  yaw  on  shot  22.  By 
repositioning  the  muzzle  0.17  m  farther  from  the  target  (shifting  the  fit  to  the  left  on  the 
graph)  for  shot  22,  the  projectile  will  be  more  likely  to  impact  the  target  during  a  portion 
of  its  yaw  cycle  where  maximum  values  occur. 
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The  length  of  the  yaw  cyele  can  be  calculated  by  taking  the  inverse  of  the  turning  rate 
and  multiplying  it  by  one  revolution.  For  example,  the  M855  fast  mode  turning  rate  is 
211  deg/m;  taking  1/211  m/deg  and  multiplying  it  by  360  deg  determines  that  the  yaw 
cycle  is  1.7  m.  Therefore,  minimum  impact  angles  can  be  expected  to  occur  for  muzzle- 
to-target  separation  distances  of  1.7  m,  3.4  m,  5.1  m,  etc.;  while  maximum  impact  angles 
can  be  expected  to  occur  for  muzzle-to-target  separation  distances  of  2.5  m,  4.2  m,  5.9  m, 
etc. 

Another  useful  total  angle  distribution  that  could  be  extracted  from  this  approach  is  the 
maximum  angle  distribution.  Using  the  linear  aerodynamics  motion  theory,  the  local 
maximum  angle  of  attack  for  an  individual  shot  is  the  sum  of  the  fast  and  slow  mode 
arms,  Kp  +  K^.  By  compiling  these  values  for  a  set  of  shots  concurrently  with  the  target 

impact  angle,  additional  characterization  of  fleet  yaw  can  be  generated.  Such  additional 
information  could  not  be  obtained  without  a  data  set  that  is  conducive  to  linear 
aerodynamics  motion  modeling. 


4.  Error  Estimation 


An  analysis  was  undertaken  to  estimate  the  accuracy  of  the  target  impact  angle,  yr , 
obtained  for  the  experiments  in  the  current  study.  Two  aspects  of  the  error  estimation 
were  considered.  First,  the  errors  of  the  trajectory  angles,  //^and^^  ,  as  determined  from 

the  linear  regression  fits  of  the  projectile  lateral  displacements  at  the  four  radiograph 
stations,  were  estimated.  Second,  the  error  of  the  target  impact  angle,  yr ,  as  determined 
from  the  linear  aerodynamics  motion  fit  of  the  pitch  and  yaw  data,  was  estimated.  In 
each  case,  a  statistical  analysis  of  the  actual  RMS  values  from  the  set  of  shots  was 
compared  to  the  distribution  of  RMS  values  from  Monte  Carlo  simulations.  The 
comparison  of  RMS  errors  allowed  a  determination  to  be  made  of  the  expected 
measurement  errors,  which  are  related  to  the  errors  of  the  quantities  of  interest. 


4,1  Trajectory  Angle  Error 

The  trajectory  angles  and 7^  ,  for  each  shot,  were  determined  via  a  linear  regression  of 
the  projectile  c.g.  lateral  displacements  as  measured  at  the  four  radiograph  stations.  The 


RMS  values  of  the  fits,  7?i ,  7?2  > 


Rp ,  are  calculated  as  7?,  =  J— 


,  where 


dj^j  is  the  displacement  of  the  projectile  at  the 7*  station,  and  is  its  lateral 

displacement  error  produced  by  the  least  squares  fit  of  the  displacement-range  pairings 
measured  at  the  four  stations.  A  frequency  histogram  was  generated  from  421  shots, 
showing  the  relative  number  of  occurrences  of  ranges  of  RMS  values  as  a  function  of 
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RMS  levels.  A  statistical  analysis  was  then  initiated  to  reproduce  this  histogram  via 
simulation  in  order  to  relate  the  RMS  histogram  to  displacement  measurement  error,  and 
subsequently  to  trajectory  angle  error. 

The  statistical  analysis  was  performed  as  follows.  It  can  be  shown  (see  appendix  A)  that 
the  maximum  likelihood  estimator  for  the  displacement  standard  deviation,  denoted  by 
|2  ^  ^ 

a ,  is  d  =  J —  YR,  .For  the  horizontal  RMS  values  from  the  actual  firings,  we  obtain 

a  standard  deviation  estimate  of  =1.31  mm;  while  in  the  vertical  direction,  we  have 
dy  =  1.21  mm. 


To  verify  that  the  model  used  to  develop  this  estimate  of  the  displacement  standard 
deviation  is  tenable,  we  compare  in  figure  14  the  relative  frequency  histogram  of  421 
horizontal  RMS  values  from  the  actual  firings  with  that  of  500,000  simulated  RMS 

values,  whose  distribution  is  Weibull'  with  a  fixed  scale  parameter  =  .9255  mm  and 

V2 

shape  parameter  2.  Figure  15  is  a  similar  comparison  for  the  vertical  RMS  values. 


0  0.5  1.0  1.5  2.0  2.5  3.0  3.5 

Horizontal  RMS  (mm) 


Figure  14.  Comparison  of  relative  frequency  histograms  of  experimental  RMS  and  simulated 
RMS  under  assumption  of  constant  errors  in  the  horizontal  direction. 


'  Several  parameterizations  exist  for  the  two-parameter  Weibull  distribution.  Compare  Mood,  Graybill  and  Boes 
(1974),  Casella  and  Berger  (1990),  and  Everitt  (1998).  In  this  paper,  we  use  the  one  advocated  by  Everitt,  that  is 

(x)  =  ^  x^  *e  >  since  this  is  the  same  parameterization  used  in  MATLAB®.  Under  this  parameterization. 


a  >  0  is  the  scale  parameter  and  h  >  0  is  the  shape  parameter. 


17 


0.16  r 

0.14 

0.12 


o  0.10 
ci> 


0.08 


ro 

aj  0.06 


0.04 


0.02 


I  Experimental 
•  Simulated  (constant 
error  assumption) 


0.5  1.0  1.5  2.0  2.5 

Vertical  RMS  (mm) 


3.0 


3.5 


Figure  15.  Comparison  of  relative  frequeney  histograms  of  experimental  RMS  and  simulated 
RMS  under  assumption  of  eonstant  errors  in  the  vertieal  direetions. 


It  is  evident  in  both  figures  14  and  15  that  the  statistical  model  used  to  estimate  the 
displacement  errors  is  not  consistent  with  the  experimental  data.  In  particular,  the  mode 
(i.e.,  the  position  of  the  peak  on  the  abscissa)  of  each  distribution  is  noticeably 
overestimated  by  the  simulated  distribution.  The  inconsistency  was  produced  by  the 
assumption  of  a  constant  standard  deviation  of  the  error  in  the  projectile  lateral  c.g. 
displacements.  A  review  of  the  complete  set  of  the  radiograph  film  shows  that  the 
images  from  different  stations  vary  noticeably  in  clarity.  Those  with  more  sharpness  and 
contrast  could  be  expected  to  produce  more  accurate  readings,  and  so  it  is  reasonable  to 
assume  that  the  displacement  error  standard  deviation  behaves  as  a  random  variable. 
Additionally,  the  processor  of  the  film  could  change  on  a  shot-to-shot  basis,  causing  the 
films  to  be  much  darker  or  lighter.  Hence,  we  incorporated  this  observation  into  the  error 
analysis  by  simulating  sets  of  RMS  values,  where  the  z*  RMS  is  drawn  from  a  Weibull 

distribution  with  varying  scale  parameter  and  shape  parameter  2. 

V2 

The  choice  of  a  distribution  to  use  in  modeling  the  displacement  errors  (that  is, 

CTj ,  0-2 , . . . ,  cr^  )  is  quite  subjective.  In  Bayesian  statistical  analyses,  the  inverted  gamma 

distribution  is  frequently  used  as  a  prior  distribution  for  variances;  under  this  precedent, 
we  considered  the  square  root  of  an  inverted  gamma  (RIG)  distribution  to  model 
displacement  errors.  We  also  considered  a  gamma  distribution,  since  this  relatively 
simple  distribution  is  quite  versatile  for  modeling  random  variables  that  are  right-skewed, 
as  would  be  expected  with  displacement  errors.  In  setting  the  parameters  of  the  gamma 
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and  RIG  distributions  -  both  have  two  parameters  -  we  imposed  the  constraint  that  the 
expected  value  equals  the  maximum  likelihood  estimate  of  the  displacement  error  under 
the  constant  error  assumption. 


The  expected  value  of  a  RIG  random  variable  with  scale  parameter  a  and  shape 


parameter  b  can  be  shown  to  equal 
b  to  vary  while  setting  a  = 


1  r(h-l/2) 
VT  r(h) 

Tib-yiy^ 


(see  appendix  B).  Hence,  we  allowed 


r(;>) 


,  and  chose  the  RIG  distribution  that  best 


matched  the  experimental  and  simulated  RMS  values.  To  objectively  select  a  set  of 
“besf  ’  RIG  parameters,  we  compared  the  empirical  cumulative  distribution  functions  of 
experimental  and  simulated,  using  the  maximum  vertical  distance  between  these  curves 
as  a  measure  of  the  goodness  of  fit,  similar  to  a  Kolmogorov-type  test  statistic  used  in 
nonparametric  distribution  fitting  (see  Conover,  1980).  Although  not  shown  here,  even 
the  “best”  RIG  distributions  had  a  much  greater  likelihood  of  high  RMS  values  than  was 
observed  in  the  experimental  results. 


The  expected  value  of  a  gamma  random  variable  with  scale  parameter  a  and  shape 
parameter  b  is  simply  their  product  ab.  Therefore,  we  again  allowed  b  to  vary  while 
setting  a  =  (j/h ,  and  chose  the  gamma  distribution  that  yielded  the  best  match  between 

the  experimental  and  simulated  RMS  values  using  the  Kolmogorov  vertical  distance  to 
evaluate  various  models. 


For  the  horizontal  data,  the  best  gamma-distributed  fit  was  obtained  with  a  scale 
parameter  of  .4027  and  a  shape  parameter  of  3.25;  for  the  vertical  data,  a  scale  parameter 
of  0.1684  and  a  shape  parameter  of  7.2  produced  the  best  fit.  Figures  16  and  17  show  the 
relative  frequency  histograms  of  the  experimental  RMS  values  overlaid  with  the  relative 
frequencies  of  values  obtained  through  simulation  with  these  “best  fit”  parameters.  The 
gamma-distributed  model  horizontal  RMS  values  still  had  relative  frequencies  that 
diminished  at  lower  RMS  values  than  observed  in  the  experimental  results,  although  not 
as  severely  as  with  the  best  RIG-distributed  model.  The  RMS  frequency  peak  locations 
agree  noticeably  better  using  the  gamma-distributed  model. 
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Figure  16.  Comparison  of  relative  frequency  histograms  of  experimental  RMS  and  simulated 
RMS  under  assumption  of  gamma-distributed  (a=.4027,  b=3.25)  errors  in  the 
horizontal  direction. 
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Figure  17.  Comparison  of  relative  frequency  histograms  of  experimental  RMS  and  simulated 
RMS  under  assumption  of  gamma-distributed  (a=.1684,  b=7.2)  errors  in  the 
vertical  direction. 
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As  we  imposed  in  choosing  the  gamma  distribution  parameters,  the  simulated 
distributions  produce  expected  standard  deviation  estimates  in  the  horizontal  and  vertical 
lateral  directions,  and  (7^  ,  of  1.21  and  1.31  mm,  respectively.  For  estimation 

purposes,  we  will  use  an  approximate  value  of  1.25  mm.  Using  this  estimate  of  the 
lateral  displacement  measurement  error,  the  Monte  Carlo  simulation  shows  the  expected 
standard  deviation  in  calculated  trajectory  angle  to  be  approximately 
0.14  deg,  as  shown  in  figure  18. 


Figure  18.  Estimated  standard  deviation  of  trajectory  as  a  function  of  reading  error. 


The  target  impact  yaw  error  analysis  to  be  presented  subsequently  demonstrates  that 
minimization  of  the  trajectory  angle  error  is  an  important  aspect  of  the  experiment  setup. 
A  review  of  the  setup  and  consideration  of  the  major  contributors  to  trajectory  angle  error 
was  performed,  and  engineering  controls  were  considered  for  future  implementation. 

The  first  consideration  is  the  manner  that  the  tube  heads  were  aligned  into  the  orthogonal 
planes.  The  heads  were  aligned  using  hand  tools,  levels  and  plumb  bobs,  and  this 
approach  provides  lateral  measurement  accuracy  to  about  2  to  3  mm.  Additionally,  in 
this  series  of  shots,  the  radiograph  fixture  was  placed  on  a  table  which  was  leveled  using 
shims.  In  an  ideal  setting,  the  fixture  would  have  the  ability  to  be  anchored  directly  to  the 
floor,  by  use  of  adjustable  legs  to  level  the  fixture.  The  final  engineering  control  measure 
would  be  the  surveying  in  of  the  tube  head  and  fiducial  wires  so  the  actual  location  of 
these  could  be  known  to  within  a  single  millimeter.  If  the  exact  locations  of  the  tube 
heads  are  recorded,  the  digitizing  program  can  adjust  the  apparent  locations  to  the  actual 
location  more  accurately.  If  these  engineering  controls  are  done  the  RMS  error  would  be 
reduced  significantly  to  only  the  placing  of  the  template  on  the  image.  In  effect,  the  film 
reading  error  would  become  a  larger  portion  of  the  total  measurement  error. 
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4,2  Target  Impact  Angle  Error 

To  estimate  the  error  of  the  target  impaet  angles,  yj,  first  the  projectile  angular  orientation 
measurement  error,  cr^ ,  needed  to  be  accurately  estimated.  Two  different  methods  of 

evaluation  were  completed  and  compared.  The  first  method  was  conducted  using  linear 
aerodynamics  Monte  Carlo  simulations  of  prescribed  pitch  and  yaw  motions  using 
prescribed  pitch  and  yaw  errors.  The  RMS  values  of  the  simulations  were  then  compared 
to  those  produced  from  the  set  of  linear  aerodynamics  fits  of  a  large  number  of  actual 
shots.  The  second  method  used  a  statistical  analysis  to  reproduce  a  histogram,  via 
simulation,  similar  to  the  RMS  values  produced  from  the  large  number  of  actual 
experimental  linear  aerodynamic  fits  to  estimate  the  error.  Upon  completion  of  the 
individual  measurement  error  assessment,  a  Monte  Carlo  simulation  was  conducted  to 
relate  the  individual  measurement  errors  to  the  target  impact  angle  error. 

The  linear  aerodynamics  Monte  Carlo  simulations  were  performed  as  follows.  Several 
specific  M855  motions  were  prescribed  by  specifying  the  six  linear  aerodynamic  motion 
parameters:  Kp  ,  ,  (/)p^  ,  (j)^^ ,  (/)'p ,  and  (j)'^  .  The  M855  turning  rates,  (/)'p  and  (j)'^  ,  were 

prescribed  to  be  their  nominal  values  for  each  simulation.  The  two  modal  arm 
magnitudes,  Kp  and  ,  were  both  fixed  at  2  deg.  The  slow  arm  reference  orientation 

was  prescribed  to  be  29  deg.  The  fast  arm  reference  orientation  was  varied  for  each 
simulation  to  position  the  specified  motion  within  a  variety  of  fast  arm  phases.  The  fast 
arm  reference  orientation  was  found  to  have  no  noticeable  effect  on  the  results  to  be 
presented. 

Next,  simulated  values  of  pitch  and  yaw  angle  were  randomly  generated  at  each  of  the 
four  radiograph  stations.  Using  the  aerodynamic  motion  parameters: 

Kp  ,  K^ ,  (/)p^  ,  (j)^^ ,  (f)'p ,  and  (j)'^  ,  as  previously  prescribed,  allowed  the  calculation  of  the 

corresponding  values  of  pitch  and  yaw  angle  at  each  of  the  four  radiograph  stations. 
Randomly  generated  values  of  pitch  and  yaw  angle  at  each  radiograph  station  were 
specified  to  be  normally  distributed  about  the  calculated  values,  i.e.,  with  no  bias.  The 
standard  deviation  error  of  the  prescribed  pitch  and  yaw  angles,  denoted  cr^ ,  were  varied 

for  three  different  cases  to  be  0.1,  0.2,  and  0.3  deg.  The  angle  error  standard  deviations 
were  assumed  to  be  the  same  in  the  pitch  and  yaw  planes.  In  each  case,  the  randomly 
generated  pitch  and  yaw  angle  data  sets  were  fitted  to  the  linear  aerodynamics  equations 
of  motions  using  the  same  least  squares  approach  used  to  analyze  data  sets  from  actual 
firings. 

For  each  prescribed  motion,  the  Monte  Carlo  simulation  was  performed  for  1000  cases  in 
which  pitch  and  yaw  angles  were  randomly  generated.  In  each  case,  the  RMS  of  the  fit, 
comprised  of  the  differences  between  the  pitch  and  yaw  angles  from  the  prescribed  and 
randomly  generated  cases,  was  recorded.  The  results  of  three  representative  simulations 
demonstrated  an  approximate  match  of  .15  deg.  Thus  the  Monte  Carlo  simulations  were 
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reevaluated  at  five  different  angle  error  standard  deviation  levels  of.l,  .13,  .15,  .18,  and 
.2  deg,  and  again  the  fit  RMS  error  was  recorded.  Figure  19  depicts  the  RMS  error  for  all 
five  Monte  Carlo  simulations  and  the  actual  experimental  RMS  errors  for  287  shots 
observed  versus  the  relative  frequency. 


RMS  Error  verses  Relative  Frequency 


.10  deg  — ♦ —  .20  deg  Experimental  Data  — • - .15  deg  — ■ —  .18  deg  — a —  .13  deg 


Figure  19.  Least-squares  fit  RMS  distribution  from  Monte  Carlo  simulations. 

The  distribution  peak  location  for  the  0.15  deg  error  simulation  appears  to  be  near  that  of 
the  actual  firings.  However,  the  relative  frequencies  of  the  0.15  deg  error  simulations 
appear  to  diminish  at  lower  RMS  values  than  the  actual  firings;  the  relative  frequencies  of 
the  0.20  deg  error  simulation  appear  to  diminish  in  a  manner  more  consistent  with  the 
actual  firings.  It  can  be  speculated  that  the  difference  between  the  distribution  shapes  of 
the  0.15  deg  error  simulation  and  the  actual  firings  is  attributable  to  the  difference  in 
magnitude  of  yaw  levels  that  each  represents.  Specifically,  the  simulations  use  a 
prescribed  maximum  yaw  of  4  deg  throughout.  The  287  actual  firings,  on  the  other  hand, 
include  approximately  60  shots  having  maximum  yaw  levels  greater  than  10  deg,  and 
several  between  18  to  22  deg.  These  relatively  high-yaw  shots  can  be  expected  to  have 
larger  RMS  fit  values,  even  though  the  ratio  of  RMS  value  to  maximum  yaw  level  may 
remain  constant  or  even  decrease  compared  to  the  low  yaw  shots.  Given  these 
considerations,  a  measurement  error  between  0.15  and  0.20  deg  is  determined  using  the 
method  of  using  linear  aerodynamics  Monte  Carlo  simulations  to  estimate  measurement 
error. 
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The  second  method  used  to  quantify  the  error  was  performed  by  statistically  fitting  the 
three-dimensional  trajectory  model  RMS  values  to  a  distribution.  This  method  lacks  the 
theoretical  basis  used  in  the  one-dimensional  trajectory  angle  case  previously  discussed, 
because  the  distribution  of  the  eight  individual  errors  which  comprise  the  RMS  follow  in 
the  complex  non-linear  equation  used  to  fit  the  data  is  unknown.  However,  the  following 
analysis  continues  to  make  use  of  the  usual  regression  assumption  that  the  sum  of  squared 
deviations  (observed  minus  predicted  values)  follow  a  chi-squared  distribution.  The 
degrees  of  freedom  associated  with  this  distribution  are  calculated  as  n-  p  =  4  ,  where 
n  =  8 ,  the  number  of  observations,  and  />  =  4  ,  the  number  of  parameters  fit  in  the  model. 

It  can  be  shown  (see  appendix  C)  that  the  maximum  likelihood  estimator  for  the  standard 

deviation  of  individual  errors  in  the  3-D  model,  denoted  by  cr ,  is  still  a  = 

Based  on  287  observations,  we  obtained  a  standard  deviation  estimate  of  <7  =  .20  deg.  A 
comparison  of  the  relative  frequency  histograms  of  287  experimental  3-D  RMS  values 
and  500,000  simulated  3-D  RMS  values  under  the  assumption  of  a  constant  standard 
deviation  of  (J  =  .20  deg  is  displayed  in  figure  20.  We  see  that  the  fit  is  poor,  and  so  we 
try  to  fit  a  model  in  which  the  standard  deviation  is  a  random  quantity  that  changes  with 
each  shot.  The  choice  of  parameters  which  minimized  the  vertical  distance  between  the 
experimental  and  simulated  RMS  cumulative  distribution  functions  was  a  =  .0164  and 
b  =  12.25 .  Simulated  RMS  values  from  this  distribution  have  an  expected  standard 
deviation  of  .20  deg.  the  same  as  that  obtained  using  the  chi-squared  distribution.  Figure 
21  compares  the  experimental  and  simulated  relative  frequency  histograms  of  the  3-D 
RMS  values  when  the  standard  deviations  are  modeled  with  this  gamma  distribution. 
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Figure  20.  Comparison  of  relative  frequency  histograms  of  experimental  RMS  and  simulated 
RMS  under  assumption  of  constant  errors. 
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Figure  21.  Comparison  of  relative  frequency  histograms  of  experimental  RMS  and  simulated 
RMS  under  assumption  of  gamma-distributed  (a=.0164,  b=12.25)  errors. 


The  linear  aerodynamies  Monte  Carlo  simulations  were  also  used  to  generate  target 
impact  angle  error  distributions,  denoted  here  as^^ ,  shown  in  figure  22.  A  statistical  test 

was  applied  to  verify  that  the  distributions  were  half-normal  distributions.  The 

estimated  standard  deviation  centered  about  zero,  denoted  cr^ ,  was  calculated  for  each 
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simulation  using  equation  13. 


cr 


r 


(13) 


For  values  of  cr^,  of  0.1,  0.2,  and  0.3  deg,  the  eorresponding  values  of  cr^  were  found  to 

be  0.096,  0.192,  and  0.286  deg,  respeetively,  and  graphieally  in  figure  23.  For  praetieal 
purposes,  the  relationship  is  linear  and  cr^  is  nearly  equal  to  cr^  for  the  experimental 

setup  of  this  study.  The  faet  that  cr^  is  nearly  equal  to  cTq,  for  the  experimental  setup  of 

this  study  is  eoineidental.  Certain  variations  in  setup  parameters  would  deerease  the 
value  of  (7^  eompared  to  cr^  .  Examples  would  be:  (1)  moving  the  target  farther 

downrange  from  the  fourth  radiograph  station,  (2)  inereasing  the  distanee  between 
radiograph  stations,  and  (3)  using  three  radiograph  stations  instead  of  four.  Multiple 
Monte  Carlo  simulations  varying  the  fast  mode  referenee  value  showed  variations  in 
(7^  to  be  negligible. 
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Figure  22.  Monte  Carlo  simulation  of  impact  yaw  error  distribution. 
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Figure  23.  Estimated  standard  deviation  of  total  impact  yaw  as  a  function  of  reading  error. 

The  linear  aerodynamies  Monte  Carlo  simulations  show  the  estimated  projectile  angle 
measurement  error  to  be  between  0.15  and  0.20  deg.  The  statistical  fit  shows  the 
measurement  error  to  be  approximately  0.20  deg.  Figure  23  shows  that  these  levels  of 
projectile  measurement  error  correspond  to  a  target  impact  angle  error  of  approximately 
0.15  to  0.20  deg. 

The  estimated  target  impact  yaw  error  includes  some  influence  of  trajectory  angle  error 
via  consideration  of  the  RMS  fit  errors  of  the  actual  firings.  However,  the  linear 
aerodynamics  model  does  not  include  this  effect,  which  is  manifested  as  a  shot-to-shot 
bias  in  the  assumed  origin  of  the  pitch  versus  yaw  coordinate  system.  That  is,  the  origin 
of  the  pitch  and  yaw  data  for  a  particular  shot  is  expected  to  be  shifted  by  the  estimated 
trajectory  angle  error,  0.14  deg.  The  effect  of  this  biased  trajectory  angle  measurement 
error  on  the  target  impact  yaw  error  can  be  conservatively  estimated  by  calculating  an 
adjusted  target  impact  yaw  error,  cr^ ,  i.e.. 


in  which  no  correlation  is  assumed  to  exist  between  cr^  and  cr^  .  For  =  0.14  deg  and 
cr^=  0.15  deg,  cr^  is  calculated  to  be  approximately  0.21  deg.  For  cr^  =  0.14  deg  and 
cr^=  0.2  deg,  cr^  is  calculated  to  be  approximately  0.24  deg. 
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5.  Conclusions 


At  the  time  of  this  writing,  approximately  500  firings  had  been  conducted  and 
subsequently  analyzed  using  the  Yaw  at  Impact  program  described  in  this  report.  The 
radiograph  measurement  technique  and  the  linear  aerodynamics  model  provide  a  cost 
effective  and  efficient  approach  to  study  the  necessarily  large  sample  sizes  required  to 
characterize  the  terminal  ballistics  of  this  class  of  small-caliber  ammunition. 

The  study  demonstrated  two  major  benefits.  First,  the  approach  provides  a  means  to 
accurately  relocate  the  gun  to  a  yaw  node  (i.e.,  minimum  yaw),  a  maximum  yaw  location, 
or  somewhere  in  between.  Second,  any  dataset  could  contribute  to  a  fleet  yaw 
characterization  because  a  reduction  for  each  shot  yields  values  of  the  fast  and  slow 
modal  arms,  Kf  and  Ks ,  whose  sum  is  the  maximum  yaw.  This  would  provide  an 
alternative  to  fleet  yaw  characterization  in  the  spark  range  after  the  initial  aerodynamics 
characterization  is  completed. 

Error  estimation  was  an  important  aspect  of  the  current  effort.  Error  levels  were 
estimated  by  relating  the  RMS  values  of  fitted  trajectory  parameters  to  those  generated 
using  statistical  simulations.  By  doing  so,  the  estimated  measurement  errors  of  the 
projectile  lateral  displacement  and  angular  orientation  data  points  could  be  obtained. 
These  error  levels  were  directly  related  to  the  trajectory  angle  and  impact  angle  errors. 

Eor  the  current  set  of  firings,  the  estimated  error  of  the  projectile  trajectory  angle  (that  is 
the  velocity  vector  orientation  relative  to  the  fiducial  wire)  was  found  to  be 
approximately  0.14  deg.  The  estimated  error  of  the  target  impact  angle  using  linear 
aerodynamics  theory  was  found  to  be  between  0.15  and  0.20  deg.  If  the  trajectory  angle 
error  is  considered  a  bias  in  the  system,  then  the  target  impact  angle  error  would  be 
between  0.21  and  0.24  deg. 

A  review  of  the  experiment  setup  was  made  and  possible  improvements  were  suggested. 
Most  significant  perhaps  is  that  the  error  estimation  methodology  provides  a  basis  for 
evaluating  alternative  radiograph  setup  parameters.  The  approach  is  a  useful  tool  for 
designing  an  experimental  setup  that  can  minimize  error  within  the  spatial  constraints 
usually  associated  with  radiograph  equipment  and  instrumentation.  At  the  very  least, 
such  considerations  could  prevent  the  use  of  experimental  parameters  that  could  create 
larger  errors  than  were  estimated  in  the  current  study. 
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Appendix  A:  Maximum  Likelihood  Estimation  of  the  Displacement 
Standard  Deviation 


To  derive  a  maximum  likelihood  estimator  for  the  displacement  standard  deviation,  a , 
from  a  set  of  RMS  values  it  is  necessary  to  know  their  probability  density  function.  It  is 


well  established  that 


is  distributed  as  a  chi-square  random  variable  with 


two  degrees  of  freedom.  Therefore  S'  =  —  T  =  —Ri  is  a  chi  random 

cr]lj=i  a 

variable  with  two  degrees  of  freedom,  having  the  probability  density  function 

/  \  ) 

(s')  =  se  for  5  >  0  .  This  is  also  the  probability  density  function  of  a  Weibuh 
random  variable  with  scale  parameter  42  and  shape  parameter  2. 

To  derive  the  probability  density  function  off?,  we  write  R  =  g{s)  =  ,  noting  that 

cr  >  0  .  Following  the  methods  described  in  Chapter  5,  Section  5  of  Mood,  Graybih,  and 
Boes  (1974),  the  probability  density  function  off?  is  given  by 


/«W=  -^g~\r)fs(g~4r)) 


_  d  2r  fir 

~  Sj  Js 
dr  a  \  a 


2  2r] 
=  —  —  e 
(j{  (7  J 


(A-1) 


for  r  >  0  .  This  is  the  probability  density  function  of  a  Weibuh  random  variable  with 
scale  parameter  and  shape  parameter  2. 

42 

Now  suppose  that  a  random  sample  of  ARMS  values  is  drawn  from  the  Weibuh 
distribution  with  scale  parameter  and  shape  parameter  2.  The  likelihood  function  is 

42 
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L[cr)  =  Y{^e 


-2y^ 

<7 


9  iV 

4"  -—'L’-i  N 

“  /T  .  1  -■ — ■ 


i=\  <y 


.IN 


cr  ,=i 


i~\ 


(A-2) 


and  the  log-likelihood  is 


in  Lip-)  =  In  4  -  IN  In  cr  -  ) 


cr  /=i  /=i 


(A-3) 


taking  the  derivative  with  respect  to  cr , 


N 


d  /  \  1  w 

- lnl(rT)  =  -2iV— -2Ii 

da 


cr  i=i 


4yr.^ 

-2N 

- h  k 


(A-4) 


Setting  this  last  expression  to  zero  and  solving  for  cr  one  gets, 


N 


N 


0  = 


-2N 


41  r/ 

i=i 


a 


41  r/ 

«  2iV  =  ^A_ 


cr  = 


a 


N  i=i 


(A-5) 


So  the  maximum  likelihood  estimator  for  the  displacement  standard  deviation,  cr ,  is 
given  by 


2  N  2 
\N^i‘ 


(A-6) 
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Appendix  B:  Expected  Value  of  the  Square  Root  of  an  Inverted  Gamma 
Random  Variable 


To  derive  the  expected  value  of  a  random  variable  that  is  the  square  root  of  an  inverted 
gamma  (RIG)  random  variable,  its  probability  density  function  is  first  needed.  First,  let 
X  be  a  gamma  distributed  random  variable  with  scale  parameter  a  and  shape  parameter  b. 
Then  the  probability  density  function  of  X  is 


(B-1) 


Then  Y  =  g[x)  =  ^]l/X  is  a  RIG  random  variable  with  density 

/r(T)=  ^g\y)fx(g\y)) 

ay 


-f— 1^— f— / 

dy\y^]  J 


=  -2t^ 


_ 1 


(B-2) 


aT[bY 


The  expected  value  of  Y  is 


E{Y)  =  \yfy{y)dy 


y-2b-ie/^y-dy. 


(B-3) 


a  1  (£>)  0 


\y-^^e/^y"dy 
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Now  let  z  =  7  ' .  Then  j  =  z  '  and  dy  =  -z  ^dz  .  Making  a  ehange  of  variable  inside 
the  integral,  we  have, 


£(r)= 


'r(6) 


2  00 

f  2/)-2 

9  00 


Jz 


a'’r(^)o 


j z^“  ^ze  Jz 


2 


2  h  7  26-3  ^  j 

J  z  —  ze  ^  “  az 


a*r(h)2 


0 


a 


2b-3 


,b-\ 


r(6); 


■ze 


dz 


(B-4) 


Notiee  that  the  expression 


within  the  integrand  is  the  probability  density 


funetion  of  a  Weibull  random  variable  with  seale  parameter  ya  and  shape  parameter  2. 


Without  loss  of  generality,  we  ean  temporarily  refer  to  this  random  variable  as  T.  This 
faeilitates  us  in  reeognizing  the  integral  as  the  {lb  -  3)^^  raw  moment  of  T.  That  is. 


£(r)= 


(B-5) 


Simplifioation  of  this  final  expression  yields 


£(r) 


b- 


T{b) 


(B-6) 
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Appendix  C:  Maximum  Likelihood  Estimation  of  the  Three 
Dimensional  Displacement  Standard  Deviation 


To  derive  a  maximum  likelihood  estimator  for  the  3-D  displaeement  standard  deviation, 
(7 ,  from  a  set  of  RMS  values  it  is  neeessary  to  know  the  probability  density  funetion  of 

these  RMS  values.  We  assume  that  =  - - 2 -  is  distributed  as  a  ehi-square 

random  variable  with  four  degrees  of  freedom.  Therefore, 


Shj  /J 

V  cr  cr  \  8  /=1  cr 


(C-1) 


is  a  ehi  random  variable  with  four  degrees  of  freedom,  having  the  probability  density 

.  s  1  3  -■ 

funetion  fs{s)  =  e  for  5  >  0 . 

If  we  let  R  =  g(s)  =  -^S ,  the  probability  density  funetion  of  R  is  given  by 


j  Vs  ^  fVs 

=  - rfs  - r 

dr  a  cr 


Vs ifVs  Y 

- — r  e 

cr  2  cr 


32  3 
=  —re 


(C-2) 


for  r  >  0  .  This  density  funetion  does  not  mateh  any  of  the  eommonly  used  distributions. 
Now  suppose  that  a  random  sample  of  ARMS  values  is  drawn.  The  likelihood  funetion 


d  ^ 

4  V,.2 


L{cr)  =  n^rfe  = 


32" 


(C-3) 
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And  the  log-likelihood  is 


Inl(o-)  =  Ain  32  -  4  Ain  o-  — ^  +  Zln(^^)-  (C-4) 


CJ  /=1  i=\ 


Taking  the  derivative  with  respect  to  <j  , 


N 


—\nL{a)  =  -AN—-A±r; 
da  a  i=i 


N 


-4A 


i=\ 


(C-5) 


Setting  this  last  expression  to  zero  and  solving  for  a  one  gets, 


N 


0  = 


-4A 


8Zt^ 


N 


/=1 


— —  4A  = 

a  a 


2  N  , 

cr=  I — Z^-^ 
\Nf:x 


(C-6) 


So  the  maximum  likelihood  estimator  for  the  displacement  standard  deviation,  a ,  is 
given  by 


2  «  2 

a  =  .  — Zf-  • 


(C-7) 
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Abbreviations  and  Symbols 


c 

pitehing  moment  eoeffleient 

7f-faetor 

Magnifieation  faetor 

Kf 

Fast  modal  arm  magnitude 

Ks 

Slow  modal  arm  magnitude 

X 

downrange  loeation  of  projeetile 

Xo 

the  designated  downrange  loeation  of  projeetile  where  x  =  0 

R 

Root  Mean  Square  (RMS) 

a 

seale  parameter  of  either  a  gamma  or  Weibull  distribution  funetion 

b 

shape  parameter  of  either  a  gamma  or  Weibull  distribution 
funetion 

Greek  Symbols 

a 

P 

Y 

T1 

aA  and  Pa 
ac  and  Pc 
ar and  Pt 
Yt 

I 

^Cen  and  Pcen 
^Fo 


pitch  angle 
yaw  angle 
total  angle 
trajeetory 

apparent  piteh  and  yaw 
eorreeted  piteh  and  yaw 
piteh  and  yaw  at  target 
total  angle  at  target 

eomplex  angle  of  attaek 
fast  arm  piteh  and  yaw  eenter 

referenee  phase  angle  of  the  fast  epieyelie  modal  arm,  evaluated  at 
x=X(t  and  real 


37 


reference  phase  angle  of  the  slow  epicyclic  modal  arm,  evaluated 
at  x=  xo  and  real 

turning  rate  of  the  fast  epicyclic  modal  arm,  real 
turning  rate  of  the  slow  epicyclic  modal  arm,  real 
standard  deviation  of  reading  error 
standard  deviation  of  error  at  target 

error  for  each  station;  difference  between  the  ac  and  Pc  angles 

obtained  from  the  radiographs  and  those  obtained  from  the  linear 
aerodynamics  free-flight  equations  of  motion 

error  on  target  caused  by  measurement  error  of  radiographs 

standard  deviation 
estimate  of  standard  deviation 
gamma  function  evaluated  at  x 
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